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Abstract 

We present results of an exact diagonalization calculation of the spectral 
function A{k, lu) for a single hole described by the t-J model propagating 
on a 32-site square cluster. The minimum energy state is found at a crystal 
momentum k = §), consistent with theory, and our measured disper- 
sion relation agrees well with that determined using the self-consistent Born 
approximation. In contrast to smaller cluster studies, our spectra show no 
evidence of string resonances. We also make a qualitative comparison of the 
variation of the spectral weight in various regions of the first Brillouin zone 
with recent ARPES data. 



PACS: 71.27.+a, 74.25.Jb, 75.10.Jm 
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The t—J model has received a lot of attention in recent years. It is believed to be the 
simplest strong-coupling model of the low energy physics of the anomalous metallic state of 
high-temperature superconductors . The Hamiltonian of the model is 

H = -t (CiaCjv + H.c.) + J £)(S< • Sj - ^TUTlj), (1) 

(ij) 

where (ij) denotes nearest neighbor sites, and c\ a , c ia are the constrained operators, c ia = 

Although understanding the t—J model doped with many holes is an important issue 
in the (potential) resolution of the high T c mystery, the single-hole state is by itself an 
interesting and important problem. For example, it leads to the many-body wave functions 
which are the starting point of any rigid-band filling analysis 0. Further, the single- 
hole model has been studied in great detail by various analytical and numerical theoretical 
techniques ||, and recently angle-resolved photoemission (ARPES) data for the insulating, 
antiferromagnetically-ordered Cu0 2 planes in Sr 2 Cu0 2 Cl 2 have become available ||, thus 
allowing for detailed comparisons between theory and experiment. To be specific, these 
results on the properties of a single hole in the Cu0 2 plane provide a direct test of how well 
the t—J model (or any other microscopic Hamiltonian) describes the low energy physics of 
the Cu0 2 plane ||. 

In this paper we report the first exact diagonalization results, found using the Lanczos 
algorithm, for a single hole described by the t—J model on a 32— site square lattice. We use 
t as the unit of energy, i.e., t — 1. Figure |l] shows the distinct k points in the reciprocal 
space of the 32-site square lattice. Previous calculations for this model were mostly done on 
the 16-site (4 x 4) square lattice, where the k points along the antiferromagnetic Brillouin 
zone (ABZ) edge (from (0, n) to (tt, 0)) are degenerate. Other square lattices that have been 
studied (18-, 20-, and 26-site) do not have the important k points along the ABZ edge, viz., 
the single-hole ground state wavevector (f,f) nor many points along the (1,1) direction 
(from (0,0) to {tt, it)). The 32-site square lattice is the smallest one which has these high 
symmetry points, and does not have the spurious degeneracy of the 4x4 square lattice. 
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Thus, this paper represents a major advance in the exact, unbiased, numerical treatment of 
an important strong-coupling Hamiltonian. 

In order for us to complete the exact diagonalization on such a large lattice, we use 
translation and one reflection symmetry to reduce the total number of basis states to about 
150 million. At k = (^r> f )> 110 reflection symmetry can be used and the total number of 
basis states is about 300 million. To study the effect of finite system sizes, we will supplement 
our results with data obtained from smaller systems: the iV = 16 (4 x 4) cluster, as well as 
a 24-site (y/l8 x \/32) cluster that includes many of the important wave vectors @. 

The electron spectral function is defined by 

A(k,a;) = £ M-^UtfWu - < + E^ 1 ), (2) 

n 

where Eq and ijjy are the ground state energy and wavefunction of the model at half filling, 
respectively, and E^ -1 and V^ -1 are the energy and wavefunction of the nth eigenstate of 
the single-hole state, respectively. A(k, uj) is calculated using a continued fraction expansion 
H with 300 iterations and an artificial broadening factor e = 0.05. We obtain A(k,uj) that 
are well converged using these quantities. 

Figure |2|(a) shows A(k, uj) at J = 0.3 from (0, 0) to (tt, 7t). At (0, 0), the spectrum has a 
quasiparticle peak at uj ~ 1.34 and a broad feature at lower energies. As k moves away from 
(0, 0) along the (1, 1) direction towards (tt, 7r), spectral weight shifts from the broad feature 
to both a higher energy quasiparticle peak and the low energy tail of the spectrum. The 
quasiparticle peak increases in intensity and shifts to higher energies, reaching the valence 
band maximum at (|, |). When k goes further towards (tv, tv), the quasiparticle peak moves 
to lower energies and its intensity drops significantly. Spectral weight move towards the 
central part of the spectrum again, eventually leaving only a very small quasiparticle peak 
at uj = 1.2313, and a broad low energy structure. Figure |2|(b) shows A(k, uj) at other distinct 
k. 

From Fig. [2] one sees that along the ABZ edge A(k,uj) are qualitatively similar ||. 
They have strong quasiparticle peaks which do not disperse much. The intensity of the 
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quasiparticle peak is the largest at (n, 0). As k moves from (tt, 0) to (f , f ), the intensity of 
the quasiparticle peak decreases and is the smallest at (|, |) along the ABZ edge. This can 
be made quantitative by calculating the quasiparticle weight, which is defined by 

z = \^\o^mi 
k ' 

and is proportional to the area under the quasiparticle peak. Table 0] shows the values of 
Zk. Zk is the largest at (tt, 0), and remains large along the edge of the ABZ. Outside the 
ABZ, Zk decreases very fast, especially along the (1,1) direction, (f , f ) is a saddle point: 
Zk is a maximum along the (1, 1) direction, but a minimum along the direction from (0, n) 
to (tt,0). 

From our spectral function results we have extracted a quasiparticle dispersion relation, 
and in Fig. ^| we display the resulting E(k), which is the location of the quasiparticle peak 
of A(k,u). The numerical values of E(\t) are tabulated in Table |I[ In agreement with 
earlier (and smaller system size) results, we find that the band maximum and minimum are 
at (| , |) and (tt, 7t), respectively [fT0|| . Further, the valence band maximum is located at 
k = (§, f ), consistent with the large theoretical effort that was applied to the single-hole 
problem [fjjj. The solid line in Fig. |3| is the self-consistent Born approximation (SCBA) result 



on a 16 x 16 square lattice fly]]. The agreement between our 32-site dispersion relation and 
that of the SCBA calculation is encouraging. 

To study the finite-size effects in our results, we can compare our cluster to smaller 
system size results. Unfortunately, a comparison of the spectral function is limited by the 
availability of the particular k points in smaller lattices. When comparing A(h, uo) of the 
16- 18-, 20-, 26- [0, and 32-site lattices at (0,0) and (tt.tt) 0, we find that the 

high energy features (including the quasiparticle peak) are not very sensitive to the system 
size. However, the low energy features are smeared out and broadened in larger systems - 



this point was also made in Ref . ||13| . We also calculated the spectral function for the 24-site 
lattice J7|. Although it is not square, this cluster has the same five k points along the (1, 1) 
direction as the 32-site lattice. A(k, uj) at these points are qualitatively similar to those of 



the 32-site results, except for detailed values of the intensity. We conclude that the 24-site 
lattice is large enough to capture the essential shape of the spectral function along the (1,1) 
direction, while the 16-site lattice is too small, especially when the lower energy features are 
concerned. In particular, the well-defined secondary peaks found at (f,|) on the 16-site 



lattice, which were interpreted to be related to the "string picture" 115], are not found in 



our 32-site system. Hence we find no evidence supporting the string picture (at least for 
this value of J). Further, the single-hole energy, defined as Eh = Eq' 1 — E$ = — £"(| , |), 
is calculated at 0.1 < J < 0.8 for the 32-site system. Fitting to the form Eh — J = a + bJ v 



gives a = —3.24, b = 2.65, and v = 0.72. This is consistent with the 16-site results |L2 
a = —3.17, b = 2.83, and v = 0.73, and also with the large cluster estimate of the SCBA 
calculations | i~lf . However, v is not | as suggested in the string picture [|15[ 



Figures £| and || show the bandwidth W = £7(|-, |) — E(ii,7r) and the quasiparticle weight 
^(1,2) at different J for iV = 16, 24 and 32. Fitting to the functional forms a + bJ u is more 
difficult. The best estimates of the coefficients are shown in the graphs. (The negative y- 
intercept of W versus J is non-physical and has been accounted for by higher order processes 
in t |16| .) For J < 0.4, W can be fitted reasonably well by a linear relation in J. Also, 
we find that Z(^k^ does not scale monotonically with iV in this range of J. Consequently, 
we feel that a precise extrapolation to large iV based on available numerics is not justified. 
However, our results support the hypothesis that at the J value of physical interest (J ~ 0.3) 
^) remains non-zero for macroscopic N. 

Finally, our results also allow for a comparison to the experimental dispersion relation 
data of Wells et al. 0. It is clear that the flat dispersion relation that we find along the 
ABZ does not agree with the measured E(k), in direct contradiction to a recent theory of 
the ARPES spectra based on the t-J model |T7|]. The behavior of along the ABZ edge 



also differs from experimental results: Z^ is a minimum at k = (§, f ) along this direction, 
while the intensity of the ARPES peak is a maximum. Our exact, unbiased numerical results 
clearly demonstrate that the t-J model cannot explain the data. However, as was shown 
in previous work it is now known that a t' is required to describe hole motion via 



a single-band model in a Cu02 plane. Thus, we defer a more quantitative comparison of 
theory and the dispersion relation obtained by experiment until a future publication wherein 
the effect of t' will be presented. A qualitative comparison to the variation of spectral weight 
as a function of k found numerically to that found in the ARPES data shows one surprising 
agreement. To be specific, note that the t—J model is a strong-coupling model which when 
undoped displays long-ranged antiferromagnetic broken symmetry. Thus one expects that 
for only one hole, there should be an equivalence of A(h, uS) under a shift of a reciprocal 
lattice vector of the magnetic lattice. As seen in Fig. ^| our results clearly do not display 
this feature. However, neither do the ARPES experiments H, and thus, unlike the original 
conjecture of Ref. ||, the strength of the on-site correlations cannot necessarily be resolved 
via photoemission. 

Summarizing, we have presented the spectral function of a single hole propagating on a 
32-site square lattice described by the t—J model. These exact, unbiased, numerical data now 
serve as an acid test of analytical theories of this important strong-coupling Hamiltonian. 
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FIGURES 

FIG. 1. (a) The 32-site square lattice employed in this study, (b) The distinct k points in the 
reciprocal space of this lattice. 



FIG. 2. The electron spectral function A(k, iv) for J = 0.3. 

FIG. 3. The quasiparticle dispersion relation of the 32-site t-J model at J = 0.3. The solid 
line is from the SCBA calculation of Ref. pT| , but with a vertical offset to make i?(§, f ) agree 
with ours. 

FIG. 4. The bandwidth of the t-J model at various J for N = 16, 24, and 32. The lines 
represent the best fits with the indicated functional form. 

FIG. 5. The same as in Fig. || but for the quasiparticle weight Z^^y 
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TABLES 




TABLE I. Quasiparticle 


weight Zk and energy E(h) of the 32-site t-J model at J 


= 0.3. The 


half-fill ground state energy 


Eq is -11.329720. 




k 




EQl) 


(0,0) 


0.145200 


1.339454 




0.234239 


1.563812 


(IT 7T\ 
\ 2 ' 2> 


0.311065 


1.832213 


/ 3-7T 3n\ 
U»4/ 


0.098759 


1.509999 


(7T,7r) 


0.005499 


1.231307 


M) 


0.211046 


1.672890 


(7T.0) 


0.341747 


1.710318 


(1,0) 


0.295290 


1.619414 


/ 3tt TV \ 
V 4 '4/ 


0.320631 


1.777018 
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Leung and Gooding Fig. 1 




Leung and Gooding Fig. 2 




Leung and Gooding Fig. 3 




Leung and Gooding Fig. 
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